AFRL-RB-WP-TR-2008-3018 


HIGH-SPEED  MAGNETOHYDRODYNAMIC  FLOW 
CONTROL  ANALYSES  WITH  3-D  SIMULATIONS 


Datta  Gaitonde 

Computational  Sciences  Branch 
Aeronautical  Sciences  Division 


JANUARY  2008 
Final  Report 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
AIR  VEHICLES  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7542 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for 
any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or  corporation; 
or  convey  any  rights  or  pennission  to  manufacture,  use,  or  sell  any  patented  invention  that 
may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Wright-Patterson 
Air  Force  Base  (AFRL/WPAFB)  Public  Affairs  Office  and  is  available  to  the  general  public, 
including  foreign  nationals.  Copies  may  be  obtained  from  the  Defense  Technical  Infonnation 
Center  (DTIC)  (http://www.dtic.mil). 

AFRL-RB-WP-TR-2008-3018  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR  PUBLICATION  IN 
ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION  STATEMENT. 


*//Signature// 

DATTA  GAITONDE 
Computational  Sciences  Branch 
Aeronautical  Sciences  Division 


//Signature// 

REID  MELVILLE,  Chief 
Computational  Sciences  Branch 
Aeronautical  Sciences  Division 


//Signature// 

MATTHEW  BURKINSHAW,  Technical  Advisor 
Computational  Sciences  Branch 
Aeronautical  Sciences  Division 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 


*Disseminated  copies  will  show  “//Signature//”  stamped  or  typed  above  the  signature  blocks. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 

Final 


1.  REPORT  DATE  (DD-MM-YY) 

January  2008 


4.  TITLE  AND  SUBTITLE 

HIGH-SPEED  MAGNETOHYDRODYNAMIC  FLOW  CONTROL 
ANALYSES  WITH  3-D  SIMULATIONS 


6.  AUTHOR(S) 

Datta  Gaitonde 


3.  DATES  COVERED  (From  -  To) 

01  December  1997-31  December  2007 


5a.  CONTRACT  NUMBER 

In-house 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

0601102 

5d.  PROJECT  NUMBER 

A03S _ 

5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

0A 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Computational  Sciences  Branch  (AFRL/RBAC)  REPORT  NUMBER 

Aeronautical  Sciences  Division  AFRL-RB-WP-TR-2008-3018 

Air  Force  Research  Laboratory,  Air  Vehicles  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7542 

Air  Force  Materiel  Command,  United  States  Air  Force 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING/MONITORING  AGENCY 

ACRONYM(S) 

Air  Force  Research  Laboratory  AFRL/RBAC 

Air  Vehicles  Directorate  n.  sponsoring/monitoring  agency 

Wright-Patterson  Air  Force  Base,  OH  45433-7542  REPORT  NUMBER(S) 

Air  Force  Materiel  Command  AFRL-RB-WP-TR-2008-3018 

United  States  Air  Force 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

PAO  Case  Number:  WS  08-0072,  15  Jan  2008.  Report  contains  color. 


14.  ABSTRACT 

Magnetohydrodynamic  studies  of  high-speed  flow  control  are  described  with  emphasis  on  understanding  fluid  response 
to  specific  plasma-based  perturbations.  Detailed  analysis  is  presented  of  the  effect  of  magnitudes  and  gradients  of 
magnetic  and  electric  fields,  their  orientation  relative  to  the  velocity  vector,  ionized  region  location  and  extent,  and 
various  nondimensional  parameters.  The  balance  between  ponderomotive  force  and  heating  is  a  major  determinant  of  the 
effectiveness  through  competition  between  work  and  ohmic  dissipation,  while  viscous/inviscid  interactions  play  a  crucial 
role  by  distorting  the  velocity  field.  The  interaction  with  an  external  circuit  through  electrodes  is  relatively  efficient 
when  fluid  is  slowed  and  energy  is  extracted,  but  yields  high  boundary  layer  heating  and  loss  of  control  performance 
when  fluid  is  accelerated.  These  observations  are  employed  to  unify  results  focused  on  a  broad  range  of  objectives 
including  heat  transfer  reduction,  3-D  separation  suppression  through  momentum  transfer,  inviscid  instability  growth 
rate  modulation  and  energy  management  in  multiple  simulated  tip-to-tail  scramjet  designs 

15.  SUBJECT  TERMS 

high-fidelity  simulation,  multiphysics  analysis,  thermochemical  nonequilibrium 


16.  SECURITY  CLASSIFICATION  OF: 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

17.  LIMITATION 
OF  ABSTRACT: 

SAR 


18.  NUMBER 
OF  PAGES 

58 


19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 
Datta  Gaitonde 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 
N/A 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


Table  of  Contents 


Section  Page 

List  of  Figures .  iv 

Acknowledgements .  v 

1  Introduction .  1 

2  Governing  Equations  and  Numerical  Model .  3 

3  Heat  transfer  mitigation  -  Type  IV  Interactions .  6 

4  Separation  Control  -  Crossing  Shock  Interactions  .  12 

5  Control  of  inviscid  disturbance  growth  -  Entropy  layer .  20 

6  Energy  management  .  28 

7  Conclusion .  41 

8  References .  44 

List  of  Acronyms,  Abbreviations  and  Symbols .  47 


iii 


List  of  Figures 


Figures  Page 

1  Features  of  Type  IV  Interaction .  7 

2  Surface  Loads  Obtained  in  Type  IV  Interaction .  8 

3  Configurations  of  Imposed  Control  .  9 

4  Key  Elements  of  Features  Obtained  with  the  Two  Conductivity  Patterns  .  .  11 

5  The  Double-Fin  Configuration  and  Surface  Oil  Flow .  13 

6  Structure  of  Flow  Past  Double-Fin  Configuration .  15 

7  Schematic  of  ECBMT  .  16 

8  Surface  Oil  Flow  Patterns  in  Double-Fin  Interaction .  17 

9  Impact  of  Control  on  Surface  Flow  Parameters  in  Double-Fin  Interaction  .  .  19 

10  Inviscid  Entropy  Layer  Instability  Problem  and  Results  Without  Control  .  .  21 

11  Plasma  Environment  Imposed  to  Control  Inviscid  Entropy  Layer  Instability  .  22 

12  Effect  of  Control  on  Energy  Norms .  24 

13  Disturbance  Vorticity  Growth  with  and  Without  Control .  25 

14  Effect  of  Profile  and  Heating  on  Disturbance  Growth .  27 

15  MGD  Energy  Bypass  Configurations  Examined  .  29 

16  Plasma  and  Electromagnetic  Environment .  30 

17  Pressure  Contours  in  Axisymmetric  Device .  33 

18  Streamlines  for  Different  Cases  Examined .  34 

19  Mean  Profiles  through  Axisymmetric  Propulsion  Device .  35 

20  Centerline  Profiles  for  Various  Cases  .  37 

21  Current  Pattern  Established  in  Generator  and  Accelerator .  38 

22  Ponderomotive  Force  Vectors  .  39 

23  Contours  of  Joule  Heating  Terms  in  the  Energy  Equation .  41 

24  Balance  Between  Energy  Terms  Along  Centerline .  42 

List  of  Tables 

1  Cases  Computed.  Q=hiteraction  parameter,  /iV  =  I  Ial!  parameter  .  31 


IV 


Acknowledgements 


This  work  was  supported  by  grants  from  the  Air  Force  Office  of  Scientific  Research  and 
high  performance  computing  time  from  the  Department  of  Defense  High  Performance 
Computing  Modernization  Program  through  centers  at  the  Aeronautical  Systems  Center, 
Naval  Oceanographic  Office,  Engineer  Research  and  Development  Center  and  the  Army 
High  Performance  Research  Center. 


v 


1  Introduction 


Plasma-based  methods  may  in  principle  accomplish  flow  control  objectives  at  a  wide  range 
of  speeds.  The  physical  mechanisms  that  account  for  the  primary  interaction  depend  on 
the  details  of  both  the  configuration  of  the  electromagnetic  field  and  the  properties  of  the 
flow.  Plasma-based  methods  are  particularly  attractive  at  high-speeds,  where  advantages 
accruing  from  the  absence  of  moving  parts,  rapid  on-demand  deployment  and  relatively 
lower  ionization  energy  budgets  are  crucial. 

In  the  low-speed  regime,  much  emphasis  has  recently  been  placed  on  radio- frequency  based 
dielectric  barrier  discharges1,2  because  of  their  ability  to  inhibit  stall  through  significant 
near  wall  electrohydro  dynamic  effects  even  at  atmospheric  pressures.  In  high-speed 
supersonic  and  hypersonic  regimes,  power  requirements  are  generally  considerably  higher 
and  the  main  techniques  typically  exploit  tailored  heat  deposition  with  microwave  or  laser 
fields,3,4  plasma  jet  injection5  and,  ponderomotive  forces  from  magnetic  fields.  This  last 
approach  is  the  focus  of  the  present  report,  and  is  denoted  magnetohydrodynamics  (MHD), 
or,  magnetogasdynamics  (MGD),  the  latter  to  reemphasize  the  fact  that  in  some  ways, 
supersonic  (kinetic  energy  dominated)  and  subsonic  or  incompressible  ionized  flows 
(pressure  dominated)  respond  in  opposite  fashion  to  imposed  magnetic  fields.6 

The  interaction  between  the  magnetic  field  and  the  fluid  depends  on  the  relative 
orientation,  magnitude  and  gradient  of  several  vector  and  scalar  fields,  including  velocity, 
U,  magnetic  induction,  B ,  current  ,j,  electric  field,  E  and  conductivity,  er,  distributions.  In 
contrast  to  control  with  microwaves  and  lasers,  which  emphasize  aerodynamic  shaping 
through  heat  deposition,  MGD  effects  include  both  force  and  energy  components.  The 
former  is  the  usual  ponderomotive  j  x  B  force,  while  energy  interactions  may  be  divided 
into  two  parts,  the  work  done  by  this  force,  and  Ohmic  dissipation  (or  Joule  heating).6,7 

Accurate  simulation  of  these  phenomena  require  integration  of  the  Navier-Stokes  and 
Maxwell  equations.  Several  recent  efforts  have  demonstrated  success  in  this  endeavor  (see 
e.g.  Refs8-11).  Typically,  a  minimal  degree  of  realism  requires  both  3-D  and  viscous  effects 
to  be  accounted  for.  However,  a  self-consistent  procedure  to  evaluate  transport  properties, 
particularly  electrical  conductivity,  is  beyond  current  capability  because  of  computational 
resources  required  and  gaps  in  knowledge  about  physical  mechanisms  that  yield  ionization. 
As  a  compromise,  the  present  calculations  employ  first-principles  for  the  fluid  dynamics, 
coupled  with  judiciously  chosen  phenomenological  approaches  for  the  electromagnetic 
environment.  Such  a  decoupling  is  reasonable  since  the  ionization  mechanism  is  anticipated 
to  be  of  the  nonequilibrium  rather  than  of  the  thermal  type.  Another  simplification  detailed 
later  is  the  use  of  the  source  term  formulation,  justified  because  the  magnetic  Reynolds 
number  is  small.  Various  theoretical  and  numerical  aspects  are  summarized  in  Section  2. 

Several  MGD-based  local  or  semi-local  control  objectives  have  been  postulated,  including 
enhancement  or  mitigation  of  i)  heat  transfer,  ii)  separation  and  iii)  evolution  of 
instabilities.  On  a  larger  scale,  magnetic  fields  have  also  been  proposed  to  reduce  drag  and 
separately  to  aid  energy  management  in  high-speed  air-breathing  propulsion  as 
incorporated  into  the  so-called  MGD  energy  bypass  strategy  (see  e.g.,  Ref.12  for  a  detailed 
discussion) . 
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In  this  article,  numerous  simulations  aimed  at  accomplishing  the  above  objectives  are  set  in 
a  unifed  context,  to  extract  generalized  insight  into  the  strengths  and  limitations  of  such 
techniques.  Localized  effects  examined  include  Type  IV  shock-on-cowl-lip  interactions  to 
reduce  heat  loads  (Section3),  crossing-shock  interactions  to  suppress  separation  (Section  4), 
and  entropy  layers  to  control  unstable  inviscid  modes  (Section  5).  On  a  larger  scale,  system 
level  objectives  of  energy  management  for  supersonic  propulsion  devices  are  explored  in 
Section  6  with  emphasis  on  an  rectangular  and  axi-symmetric  configurations.  Of  the 
distinct  strategies  proposed  in  the  literature,  here  a  generator  is  mounted  in  the  inlet  of  the 
device,  while  an  accelerator  is  located  downstream  of  the  combustor. 
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2  Governing  Equations  and  Numerical  Model 


The  aerospace  environment  of  interest  is  characterized  by  relatively  low  conductivity, 
measured  through  the  magnetic  Reynolds  number,  Ra  =  (imaUL  «  1,  where  pm  is  the 
magnetic  permeability,  U  and  L  denote  velocity  and  length,  respectively,  and  subscript  ref 
denotes  a  reference  value.  Under  this  condition,  it  is  reasonable  to  assume  that  the 
magnetic  field  induced  by  the  current  is  much  smaller  than  that  imposed  on  the  flow.6  A 
considerable  simplification  can  then  be  realized  by  using  the  low  magnetic  Reynolds 
number,  or  source  term,  formulation  of  the  governing  equations.  Details  are  presented  in 
the  literature  (e.g.,  Ref.13),  and  are  not  reproduced  here. 


The  transport  properties  are  the  molecular  viscosity  p,  obtained  with  Sutherland’s  law, 
and  the  electrical  conductivity  a.  In  addition  to  the  Reynolds  number, 

Re  =  prefUrefLref / firef,  Mach  number,  M,  and  Prandtl  number,  Pr  =  prefCp/kref,  an 
interaction  parameter  Q  =  arefB^ejLref / (prefUref)  provides  a  scale  factor  between 
ponderomotive  and  inertial  forces. 

The  electrical  quantities  specifying  the  plasma  environment  are  cr,  B ,  j  and  E.  As  noted 
earlier,  a  fully  self-consistent  approach  to  derive  these  quantities  remains  unfeasible  at  this 
time,  especially  in  three-dimensional  situations.  Consequently,  the  transport  property  a  is 
assumed  to  be  specified,  and  since  it  is  case  dependent,  is  described  in  the  context  of  each 
problem.  The  parameters  chosen  for  all  simulations  correspond  to  wind-tunnel  conditions, 
since  in  several  cases  results  are  available  (without  MGD)  to  provide  a  measure  of 
validation.  An  added  benefit  to  considering  cold  flows  is  that  the  frozen  gas  assumption 
neglecting  high-temperature  effects  is  also  then  reasonable. 


Since  the  induced  magnetic  field  is  negligible  in  the  low  Ra  approximation,  B  is  the  known 
(imposed)  induction  field.  The  current  j  is  obtained  with  the  phenomenological  form  of  the 
generalized  Ohm’s  law:14 


j  =  a  ■ 


E  +  U  x  B 


(1) 


Expressions  for  the  3-D  conductivity  tensor,  cr,  with  Hall-current  and  ion-slip  have  been 
presented  in  Ref.15  The  impact  of  these,  particularly  the  Hall  component  of  the  current, 
can  be  significant16  for  the  environment  under  consideration.  Techniques  to  minimize  the 
Hall  effect  include  the  use  of  segmented  instead  of  continuous  electrodes  as  utilized  for 
energy  management  simulations  (Section  6).  The  electric  field  E  is  determined  from  the 
current  continuity  condition: 

V'J=0.  (2) 

Introducing  a  scalar  potential,  E  =  —  \/  cf>,  the  equation  solved  is  as  follows: 


v- 


d  •  [v0] 


V- 


d-U  xB 


(3) 


Some  of  the  simulations  discussed  below  are  turbulent.  Although  the  effect  of  the  magnetic 
field  on  turbulence  is  complex,  for  the  present  purpose,  it  is  sufficient  to  consider  an 
engineering  approach  to  mimic  the  results  of  fine-scale  turbulence  on  the  mean  flow.  The 
method  utilizes  the  popular  two-equation  k  —  e  model,  where  k  is  the  turbulence  kinetic 
energy  and  e  is  its  dissipation.  Thus,  the  mean  flow  equations  are  assumed  to  be  of  the 
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same  form  as  in  the  laminar  situation,  with  the  molecular  viscosity  p  being  replaced  by  the 
sum,  n  +  fit,  where  pt  is  the  eddy  viscosity  and  a  turbulent  Prandtl  number  (. Prt  =  0.9)  is 
introduced  in  the  standard  fashion.  The  baseline  two-equation  k  —  e  model,  including 
low-Reynolds  number  terms,  employs  the  work  of  several  authors  and  has  been  described  in 
Ref.15  New  terms  are  added  to  the  original  model  to  reproduce  some  of  the  anticipated 
effects  of  the  magnetic  held  in  a  simple  yet  effective  manner  based  on  the  recent  work  of 
Kenjeres  et  al.17 

Algebraic  manipulation  permits  the  fluid  dynamic  equations,  the  potential  equation  and 
the  turbulence  equations  to  be  written  in  flux  form  as  follows: 


dX  OF,  dGj  <9iP/  ()FV  [)GV  dHv  A 

ir  +  iy  +  iy  +  iy“  +  ^r  +  ^r+s 


where  a  general  curvilinear  coordinate  transformation  has  been  introduced,  x  =  x(£,r),  £), 
V  —  V  (£,  V,  C)  and  z  =  z(£,  y,  C),  in  order  to  facilitate  the  treatment  of  complex 
configurations.  For  the  fluid  equations,  X,  is  the  vector  1/J{p,  pu,  pv,  pw,  pe},  where  J  is 
the  Jacobian  of  the  coordinate  transformation.  Fj,  Gj,  and  Hj  contain  terms  relevant  to 
inviscid  fluxes,  Fy,  Gy,  and  Hy  include  effects  due  to  viscosity,  and  S  is  the  source  term 
containing  electromagnetic  interaction  terms.  For  the  Poisson  equation,  X  is  a  scalar 
(=  0/J),  and  the  time  derivative  is  added  to  facilitate  pseudo  time  advancement  to  steady 


state.  Terms  with  subscript  V  are  zero,  and  S  is  essentially  y  •  a  (-U 


x  B 


For  the 


turbulence  equations,  AJ;e  is  the  vector  pk,pe },  and  further  particulars  of  each  term  in 
Equation  4  may  be  found  in  Refs. 10,13,18 


Since  the  influence  of  the  magnetic  held  is  restricted  to  the  source  term  in  the  present  low 
Ra  approach,  conventional  computational  fluid  dynamics  (CFD)  techniques  have  been 
incorporated.10,13  For  problems  with  strong  shocks,  the  Roe  upwind-biased  scheme19  is 
employed  to  nominal  third-order  accuracy,  together  with  a  harmonic  limiter  to  enforce 
monotonicity. 10  When  small  disturbances  are  to  be  resolved,  as  in  the  analysis  of  unstable 
modes,  the  fourth-order  compact  difference  scheme  is  utilized  together  with  a  Pade-type 
hlter  of  eighth-order.20  These  high-order  schemes  are  also  exclusively  utilized  for  spatial 
discretization  of  the  Poisson  equation. 


Two  different  time-integration  procedures  are  considered.  For  steady  state  problems,  or 
where  solid  walls  and  attendant  fine  meshes  are  encountered,  an  approximately  factored 
scheme  is  utilized21  in  diagonalized  form.22  The  method  is  extended  to  the  Poisson 
equation  as  in  Ref.23  For  unsteady  problems,  the  fluid  equations  are  solved  with  the 
classical  fourth-order  Runge-Kutta  scheme.24 


The  code  employed  to  solve  the  above  equations  has  been  verified  for  several  canonical 
configurations.  Results  for  a  flat  plate  boundary  layer  subject  to  a  magnetic  held  have  been 
successfully  compared  with  similarity  theory  in  Ref.10  Aspects  of  channel  hows  bounded  by 
segmented  electrodes  have  also  been  compared  to  theory  and  previous  simulations  by  other 
researchers  in  Ref.13  The  ability  of  the  method  to  accurately  capture  the  sharp  gradients  of 
electric  potential  and  the  consistent  build  up  of  electric  held  at  electrode-insulator 
junctures  were  successfully  demonstrated  together  with  various  effects,  such  as  held 
reversal  near  electrodes,  which  are  also  observed  in  more  complex  devices.  Another  key 
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phenomenon  anticipated  in  high-speed  flows  is  the  combination  of  events  termed  swept 
shock/boundary  layer  interactions.  A  series  of  previous  papers  (e.g.,  Refs.25,26)  have 
demonstrated  the  capability  of  the  scheme,  through  exhaustive  comparison  with 
experiment,  to  capture  the  key  features  of  such  flows,  including  3-D  separation,  attachment 
and  vortex-shock  interactions. 

The  use  of  curvilinear  coordinate  transformations  facilitates  the  refinement  of  the  mesh  in 
regions  of  high  gradients  for  superior  accuracy.  These  simulations  have  been  subjected  to 
standard  tests  to  demonstrate  a  satisfactory  degree  of  mesh  independence.  Further 
information  is  provided  in  the  context  of  each  problem. 

Boundary  conditions  employed  are  also  relatively  straightforward.  The  fluid  dynamic 
variables,  U,  p  and  T  are  treated  in  the  same  manner  as  in  the  no-MGD  situation.  On  all 
solid  surfaces,  the  no-slip  condition  is  enforced  and  the  pressure  gradient  is  set  to  zero, 
since  the  magnetic  field  gradient  is  small.  The  wall  temperature  is  specified  to  represent 
either  a  cooled  or  adiabatic  wall,  depending  on  the  problem.  The  downstream  boundaries 
are  predominantly  supersonic  and,  thus,  the  zero  gradient  condition  is  applied.  At  inflow 
boundaries,  the  flow  vector  is  specified.  Current  and  electric  field  conditions  are  established 
through  the  electric  potential  on  the  boundaries.  Details  of  this  procedure  have  been 
described  in  Ref.13  Briefly,  the  potential  at  electrodes  is  enforced  based  on  appropriately 
chosen  load  factors  (the  external  circuit  is  not  presently  solved).  At  insulators,  the  normal 
component  of  current  is  set  to  zero.  Numerical  implementation  of  this  condition  in  general 
curvilinear  coordinates  is  not  trivial  but  high-order  stable  boundary  implementations 
successfully  developed  in  Ref.13  have  been  employed. 
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3  Heat  transfer  mitigation  —  Type  IV  Interactions 


The  shock-on-shock  interaction  occurring  near  a  surface  is  examined  as  a  prototypical 
instance  where  high  heat  loads  occur.  Of  the  many  patterns  identified  by  Edney,27  the 
Type  IV  interaction  described  below  is  typically  the  harshest  because  of  the  generation  of  a 
supersonic  jet  which  impinges  on  the  surface,  often  with  catastrophic  results.  The  flow 
parameters  are  chosen  from  the  experimental  configuration  of  Ref.28  The  cowl  lip  is 
modeled  as  a  cylinder  of  radius,  R  =  0.038m,  the  Mach  number  is  8.03,  =  111.6 K, 

Ux  =  1700 m/s,  and  the  Reynolds  number,  Re  =  1.3  x  105.  The  impinging  shock 
corresponds  to  a  flow  deflection  angle  of  12.5°.  Thus  M2  =  5.25,  p2/Poo  =  7.13  and 
P2/P00  =  3.33  where  subscript  2  refers  to  conditions  downstream  of  the  shock.  Figure  la 
shows  the  grid  employed  for  the  problem.  Points  are  clustered  near  the  body  to  adequately 
resolve  the  high  near  wall  gradients. 

The  computed  flowfield  obtained  in  the  Type  IV  interaction  is  shown  in  Figure  lb.  The 
sequence  of  events  is  initiated  when  the  impinging  shock  intersects  the  cowl  shock,  a 
situation  that  arises  as  an  unintended  consequence  when  external  compression  ramp  shocks 
are  focused  on  cowl  lips  to  optimize  mass  capture.  Two  triple  points  are  formed,  TP1  and 
TP 2,  from  which  emanate  two  shear  layers.  The  supersonic  jet  arising  between  these  layers 
passes  through  embedded  shocks  and  expansions  and  finally  through  a  terminating  jet 
shock  to  impact  the  surface,  yielding  sharp  pressure  and  heat  transfer  peaks  shown, 
respectively,  as  the  curve  Baseline  in  Figures  2a  and  b,  in  which  the  re-axis  is  the  angle 
along  the  cylinder  in  degrees,  while  the  ordinate  is  the  appropriate  surface  load  normalized 
by  the  peak  value  (without  control).  Away  from  the  location  of  this  peak,  the  heat  load 
reduces  rapidly  as  the  stagnated  fluid  accelerates  around  the  body  to  supersonic  speeds, 
encountering  flow  features  which  yield  smaller  undulations.  Although  there  is  an  inherent 
unsteadiness  in  the  interaction  (see  for  example,  Ref.29),  this  analysis  is  focused  on  the 
mean  flow  -  under  most  conditions,  the  effect  of  MGD  is  in  fact  to  stabilize  the 
unsteadiness  of  the  interaction. 

For  control,  several  electromagnetic  configurations  with  and  without  electrodes  are 
considered,  only  some  of  which  are  described  here  to  illustrate  key  observations.  In  the  first 
arrangement,  Figure  3a,  no  electrodes  are  employed  and  the  magnetic  field  is 
circumferential,  decaying  in  the  radial  direction  as  is  consistent  with  an  electric  current 
along  the  axis  of  the  cylinder.  Two  different  nonequilibrium  conductivity  distributions, 
denoted  a  a  and  aB,  respectively,  are  utilized  in  separate  computations  (Figure  3a).  The 
first,  a  a,  spans  the  impinging  jet  and  adopts  the  strategy  of  slowing  down  the  high-speed 
fluid  near  the  surface.  The  second,  aBj  is  targeted  further  away  from  the  body,  and  is 
designed  to  alter  the  location  and  characteristics  of  the  primary  triple  point  TP1.  Each 
conductivity  distribution  is  generated  from  modified  Gaussians  mimicking  those  obtained 
from  e-beam  approaches,  further  details  of  which  may  be  found  in  Ref.30  An  added 
dimension  can  be  incorporated  by  mounting  electrodes  on  the  surface  of  the  cylinder  as 
shown  in  Figure  3b  in  conjunction  with  a  cr a  conductivity  distribution.  Setting  the 
electrodes  at  specific  potentials  establishes  an  imposed  electric  field  and  generates  a  current 
from  anode  to  cathode.  By  locating  the  electrodes  on  either  side  of  the  peak  load,  and 
imposing  a  spanwise  magnetic  field  (pointing  out  of  the  plane  of  the  figure),  the  j  x  B  force 
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(a)  Mesh 
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(b)  Flow  pattern  with  Mach  number 
Figure  1:  Features  of  Type  IV  Interaction 
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(a)  Surface  pressure 

D  Baseline 
-  -A-  -oA  (Qeff=1 2.8) 


(b)  Heat  transfer  rates 

Figure  2:  Surface  Loads  Obtained  in  Type  IV  Interaction 
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(a)  Electrodeless  arrangement 


Figure  3:  Configurations  of  Imposed  Control 
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can  be  assured  to  point  away  from  the  body,  thus  reducing  the  kinetic  energy  of  the  jet. 

Key  findings  on  the  control  effect  of  these  three  configurations  on  surface  pressure  and  heat 
transfer  are  shown  in  Figure  2a  and  2b,  respectively.  The  effective  interaction  parameter 
chosen  depends  on  the  sensitivity  of  the  flow  to  the  plasma  perturbation,  and  are  thus 
different  for  the  various  cases.  Close  scrutiny  of  these  figures  indicates  that  the  presence  of 
electrodes  is  detrimental  in  the  effort  to  reduce  surface  loads.  The  surface  pressure  peak 
remains  undiminished,  though  its  location  moves  nearer  the  lower  electrode  (cathode) 
location  (smaller  6).  More  striking,  the  heat  transfer  rate  shows  two  peaks,  one  at  the 
cathode,  and  another  much  larger  value  almost  twice  as  high  as  the  original  maximum  in 
the  vicinity  of  the  anode,  which  is  the  upper  electrode.  The  development  of  peaks  near  the 
electrodes  is  linked  to  the  effect  of  near  wall  Joule  heating  (j2/a).  This  dynamics  is 
dictated  by  the  requirement  that  all  current  in  the  flow  be  channeled  through  the 
electrodes,  near  which  the  current  concentration  is  largest.  Since  sheath  effects  have  been 
ignored,  the  asymmetry  in  response  near  the  anode  and  cathode  may  be  traced  to 
differences  in  induced  electric  field  arising  in  response  to  the  different  local  motional 
electromotive  force  of  the  boundary  layer  recovering  after  stagnation. 

The  electrodeless  configurations  yield  considerably  lower  Joule  heating  rates  and  the 
results  are  more  promising.  Both  electrodeless  approaches  reduce  surface  pressure  peaks  in 
a  manner  generally  proportional  to  the  interaction  parameter.  However,  the  reduction  in 
heat  transfer  rate  is  similar  with  both  approaches,  indicating  the  superiority  of  the  erg 
distribution  since  the  effective  interaction  parameter  is  smaller.  The  difference  in 
performance  may  be  connected  to  the  reaction  generated  in  the  fluid.  Since  the  walls  are 
insulated,  the  solenoidal  constraint  on  the  electric  current  causes  eddy  loops  to  appear,  as 
shown  for  the  initial  transient  in  Figure  4a  for  a  a-  This  pattern  results  from  the  velocity 
differential  across  the  shock,  yielding  a  local  current  direction  that  points  in  the  U  x  B 
direction  upstream  of  the  shock,  thus  decelerating  the  flow.  However,  the  return  path  is 
established  by  an  induced  electric  field  through  the  slower  near  wall  fluid.  Here,  flow 
accelerates  toward  the  wall,  and  the  impact  of  control  is  diminished.  With  erg  however,  the 
closed  current  field  perturbs  the  location  of  the  primary  triple  point  sufficiently  to  change 
the  observed  pattern  from  the  harsh  Type  IV  to  the  relatively  milder  Type  III  as  shown  in 
Figure  4b.  Even  though  this  configuration  requires  establishment  of  a  suitable  plasma 
environment  further  away  from  the  surface,  it  indicates  clearly  that  control  is  most  effective 
when  employed  to  alter  the  prevailing  fluid  dynamics  to  yield  a  milder  interaction. 

A  measure  of  the  physical  values  of  nondimensional  parameters  plotted  in  Figure  3  may  be 
obtained  from  the  interaction  parameter,  Q,  and  the  freestream  flow  values.  For  the 
present  case,  pref  =  /y^,  Uref  =  U^,  and  Lref  =  R ,  which  for  reasons  noted  above  are 
chosen  to  represent  wind-tunnel  rather  than  flight  conditions.  From  the  definition  of  Q, 
crrefB2ef  =  QprefUref/Lref.  If  the  reference  value  of  the  magnetic  field  is  then  chosen,  the 
reference  electrical  conductivity  can  be  extracted.  For  a  nominal  Q  —  1,  the  flow 
parameters  yield  crrefB2ef  ~  1.2 mho  ■  T2/m.  The  current  is  scaled  by  arefUrefBref.  If  a  IT 
magnetic  field  is  considered,  the  value  of  unity  for  electrical  conductivity  in  Figure  3 
corresponds  to  a  =  1.2 mho/m.  The  highest  value  of  Q  examined  requires  arefB'2ej  ~  15.4. 
Doubling  the  magnetic  field  has  the  effect  of  requiring  one  fourth  the  electrical 
conductivity  to  obtain  the  same  effect. 
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(a)  Current  pattern  with  a  a 
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(b)  Type  III  pattern  with  <jb  (Mach  contours) 

Figure  4:  Key  Elements  of  Features  Obtained  with  the  Two  Conductivity  Patterns 
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4  Separation  Control  -  Crossing  Shock  Interactions 


The  eddy  current  pattern  shown  earlier  in  Figure4a  may  be  naturally  associated  with  a 
mechanism  to  transfer  momentum  from  high-speed  regions  of  the  flow  to  lower  speed 
regions.  In  the  outer  region,  the  local  interaction  is  generator-like;  body  forces  reduce  fluid 
speed  and  work  is  done  in  establishing  the  current  (U  ■  (j  x  B)  <  0).  The  reverse  dynamic 
occurs  in  the  near  wall  low-speed  region,  where  work  is  done  in  accelerating  the  fluid 
(U  ■  (j  x  B)  >  0).  It  is  well  known  that  properly  oriented  enhancement  of  lower  speed  (near 
wall)  fluid  momentum  can  aid  in  suppressing  separation.  Thus,  MGD  can  not  only  be 
employed  to  induce  separation  (see  e.g.,  Ref.10)  but  suitable  arrangements  of  plasma 
parameters  have  the  potential  to  eliminate  separation  as  well. 

To  illustrate  this  possibility  in  a  complex  configuration,  this  eddy-current-based  momentum 
transfer  (ECBMT)  approach  is  considered  in  a  3-D  situation  where  the  relative  orientations 
between  the  velocity,  magnetic,  current,  electric  and  body  force  fields  are  not  intuitively 
obvious.  Specifically,  secondary  separation  occurring  in  the  double-fin  configuration,  shown 
in  Figure  5a,  is  explored.  The  problem  isolates  the  physics  of  instances  where  an  incoming 
equilibrium  turbulent  boundary  layer,  such  as  that  developing  on  a  fuselage,  interacts  with 
shock  waves  generated  by  two  compression  surfaces,  which  may  be  considered  to  represent 
sidewalls  of  an  inlet.  Depending  on  the  flow  parameters,  regimes  from  weak,  where  the 
incoming  boundary  layer  does  not  separate,  to  strong,  where  separation  occurs,  have  been 
characterized  in  the  literature.  The  no-control  flow  has  been  described  in  several  prior 
publications  (see,  e.g.,  Refs.25,31  and  citations  therein).  The  freestream  Mach  number  and 
stagnation  conditions  are  M  =  4.961,  PQ  =  2.2  Mpa,  and  Ta  =  427  K,  resulting  in  a 
freestream  unit  Reynolds  number  of  36.5  million  per  meter  and  the  wall  temperature  is 
Tv,  =  295  K .  These  parameters  model  the  experiments  of  Zhcltovodov,  et  al. 32,33  The  two 
fins  are  each  placed  at  23  degrees  angles  of  attack  and  symmetry  is  invoked  to  compute 
only  half  the  domain,  though  for  clarity,  some  converged  results  are  reflected  prior  to 
plotting.  Details  of  boundary  conditions  and  mesh  resolution  may  be  found  in  Ref.30 

Figure  5b  depicts  the  computed  surface  streamline  pattern.  Various  lines  of  coalescence  S 
and  divergence  R  are  observed,  some  of  which  are  marked  in  the  figure.  The  lines  (Si,Ri) 
represent  primary  separation  and  attachment.  Downstream,  the  two  primary  lines  of 
separation  evolve  from  each  side  of  the  symmetry  plane  into  a  complex  structure  yielding  a 
centerline  of  attachment,  straddled  by  two  lines  of  off-centerline  separation.  The  line  of 
coalescence,  marked  S3,  aligned  along  an  angle  approximately  midway  between  the  primary 
separation  and  attachment  lines,  indicates  secondary  separation  (SS),  which  is  the  major 
focus  of  the  present  demonstration.  There  is  also  a  clearly  identifiable  attachment  line  in 
its  close  vicinity  as  required  for  consistency.  These  lines  terminate  downstream  where  they 
form  a  critical  point  pair  (not  marked).  Although  SS  is  commonly  observed  in  laminar 
flows,  the  phenomenon  is  relatively  rare  in  shock/turbulent  boundary  layer  interactions.  In 
Ref.,25  it  has  been  shown  that  each  feature  observed  in  Figure  5b  also  arises  in  the 
experiments  of  Ref.32 

Figure  6a  depicts  the  kinematic  structure  with  judiciously  chosen  stream  ribbons.  The 
incoming  boundary  layer  separates  along  the  primary  separation  line  and  does  not  reattach, 
but  rather,  becomes  narrow  and  curved  downstream  (only  a  part  of  this  surface  is  shown  to 
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Inviscid  Shock 


(a)  Schematic  of  problem 


Figure  5:  The  Double-Fin  Configuration  and  Surface  Oil  Flow 
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visualize  the  other  regimes).  Fluid  attaching  at  the  primary  line  of  attachment,  Ri,  may  be 
distinguished  into  three  different  components.  The  vortex  interaction  regime  contains  the 
fluid  that  separates  from  the  downstream  side  of  Si.  A  vortex  filament  pair  (not  shown) 
originates  from  two  foci  on  either  side  of  the  symmetry  plane.  The  second  regime  consists 
of  fluid  that  forms  the  centerline  vortex  while  the  last  regime  comprises  flow  attaching  at 
Ri  and  separating  at  the  line  of  SS,  S3.  Fluid  attaching  at  the  centerline  originates  outside 
the  incoming  boundary  layer  from  near  the  fin  leading  edges  and  sweeps  across  the  domain 
prior  to  attachment.  The  shock  structure  is  shown  with  the  pressure  on  several  cross- flow 
planes  in  Figure  6b.  Far  from  the  plate,  the  intersection  of  the  inviscid  shock  clue  to  the  fin 
with  the  cross-flow  plane  is  clearly  observed.  The  interaction  of  the  shock  with  the 
boundary  layer  produces  the  classic  3-D  A-pattern.34  The  two  As  from  either  side  of  the 
symmetry  plane  cross  each  other  in  downstream  regions,  in  a  complex  manner  that  has 
been  described  in  Ref.26  Fluid  near  the  plate  traverses  in  a  nearly  spanwise  direction,  and 
encounters  an  adverse  pressure  gradient,  shown  later  in  plots  of  surface  pressure. 

Figure  7  shows  schematically  a  possible  approach  to  accomplish  ECBMT  for  the  double- fin 
interaction.  The  region  between  R\  and  S3  is  assumed  to  be  the  footprint  of  a  region  of 
ionization,  extending  a  few  boundary  layer  heights  above  the  plate.  The  entrainment  of 
high-speed  fluid  near  the  plate  through  the  A-structure  greatly  facilitates  the  ECBMT 
method,  which  requires  that  the  ionized  region  extend  into  high  velocity  regions.  A 
uniform  transverse  magnetic  held  is  established  and  the  plate  is  assumed  to  be  insulated. 

A  complex  fluid-plasma  interaction  is  anticipated,  depending  on  the  interplay  between  the 
velocity,  magnetic,  current,  induced  electric  and  body  force  fields.  Notional  orientations  of 
these  holds  are  marked  in  Figure  7. 

Several  computations  exploring  ECBMT  have  been  performed  and  described  in  detail  in 
Ref.30  In  each  case,  the  electrical  conductivity  is  assumed  to  be  established  in  a  region 
above  the  plate  to  different  heights.  To  highlight  the  observations,  here  only  the  case  where 
the  ionized  region  extends  5  5  above  the  plate  region  and  the  interaction  parameter  is  unity 
(denoted  H5Q1)  is  described.  Figure  8a  shows  the  simulated  surface  oil  how  with  control. 
Close  scrutiny  reveals  that  the  effect  of  MGD  is  to  eliminate  the  secondary  separation  and 
attachment  lines.  The  surface  oil  how  lines  are  more  conical  in  form,  with  an  increased 
streamwise  component  compared  to  the  no-MGD  case.  Examination  of  the  vertical  velocity 
component  above  the  plate  (not  shown  for  brevity)  indicate  that  positive  values  observed 
without  control  are  eliminated,  consistent  with  the  elimination  of  manifestations  of 
secondary  separation  in  the  surface  pattern.  Post-processing  also  indicates  that  the  eddy 
current  pattern  is  consistent  with  that  anticipated  and  depicted  in  the  schematic  of  Figure6. 
However,  further  numerical  experiments  are  performed  to  demonstrate  that  it  is  truly  the 
transfer  of  momentum  that  is  responsible  for  the  observed  elimination  of  separation. 

In  the  hrst  simulation,  to  isolate  the  impact  of  Joule  heating,  simple  heat  release  is  modeled 
on  the  same  overall  spatial  distribution  as  the  conductivity,  but  with  the  height  restricted 
to  0.5  6.  The  local  normalized  heat  release  parameter  is  set  to  20  percent  of  the  incoming 
enthalpy.  The  surface  pattern  for  this  case,  shown  in  Figure  8b,  has  clearly  worsened  from 
the  standpoint  of  separation,  compared  with  the  no-MGD  case  of  Figure  5b.  The  features 
associated  with  secondary  separation  are  now  not  only  far  more  clearly  defined,  but  also, 
new  structures  are  evident  near  the  centerline.  In  contrast,  with  MGD,  the  results  of 
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Figure  6:  Structure  of  Flow  Past  Double-Fin  Configuration 
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Figure  7a  indicate  primarily  a  local  effect  in  the  vicinity  of  the  applied  control. 

The  possibility  that  separation  is  eliminated  because  the  entrained  fluid  is  made  more 
turbulent  by  the  plasma  interaction  has  also  been  evaluated  by  examining  eddy  viscosity 
values  with  and  without  control.  The  combined  effect  of  Joule  heating  and  damping 
tendency  of  the  magnetic  held  result  in  reduction  of  eddy  viscosity  values  in  the  ionized 
region  from  0(1000)  for  the  no-MGD  case  to  essentially  negligible  levels,  indicating  a 
laminarization  of  the  local  how  held.  Clearly,  the  elimination  of  SS  is  not  associated  with 
heightened  turbulence,  but  rather  with  ponderomotive  force  effects. 

The  effect  of  control  on  surface  pressure  and  heat  transfer  is  shown  in  Figures  9a  and  b, 
respectively  along  a  cross-how  station  located  at  A"  =  40.5.  The  adverse  pressure  gradient 
inducing  secondary  separation  is  evident  at  z/8  ~  13  in  the  curve  titled  No-MGD  in 
Figure  9a.  The  case  with  MGD  shows  a  substantial  reduction  in  pressure  near  the  fin-plate 
corner  and  maximum  values  diminish  by  about  40  percent.  The  howheld  modification  also 
eliminates  the  pressure  rise  associated  with  SS  at  z/8  ~  —13,  though  as  noted  earlier,  the 
fluid  has  a  larger  streamwise  component  with  MGD  control  than  without.  Figure  9b 
exhibits  the  effect  on  heat  transfer  rates.  Since  the  quantitative  values  are  known  to  be 
significantly  overpredicted,35  only  trends  are  reported  by  normalizing  values  with  the 
no-MGD  heat  transfer  rate.  Despite  the  electrodeless  nature  of  the  MGD  arrangement,  an 
increase  is  detected  in  the  double-fin  case.  In  contrast,  the  electrodeless  a  a  approach  to 
control  the  Type  IV  interaction  (Section  3)  yielded  a  heat  transfer  reduction.  The 
difference  in  response  arises  from  configurational  details:  in  the  Type  IV  interaction,  the 
overall  effect  of  control  is  to  interfere  with  the  impingement  process  responsible  for  high 
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(a)  With  control 
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heat  transfer  rates,  whereas  in  the  present  3-D  situation,  the  principal  effect  is  acceleration 
of  fluid  parallel  to  the  surface,  after  impingement  at  the  line  of  attachment.  Not 
surprisingly,  heat  addition  by  itself  yields  higher  peaks  (Figure  9b)  than  with  MGD,  about 
35  percent  higher  than  the  no-MGD  case. 
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(a)  Cross-flow  pressure 


(b)  Cross-flow  heat  transfer 

Figure  9:  Impact  of  Control  on  Surface  Flow  Parameters  in  Double-Fin  Interaction 
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5  Control  of  inviscid  disturbance  growth  -  Entropy 
layer 

To  explore  the  effect  of  electromagnetic  control  on  the  growth  of  inviscid  instabilities,  an 
entropy  layer  with  a  generalized  inflection  point  is  employed.  The  specific  situation  occurs 
above  a  flat  plate  with  a  blunted  leading  edge  placed  in  a  Mach  6  freestream,  as  shown  in 
Figure  10a.  The  nonlinear  Euler  equations  are  solved  with  a  fourth  order 
compact-difference  and  eighth  order  Pade  filter  for  spatial  discretization,  and  a  fourth 
order  Runge-Kutta  method  in  time.  The  normalized  base  flow  properties  of  the  entropy 
layer  are  obtained  from  Ref.36  using  the  theory  of  Refs.37,38  and  are  shown  in  Figure  10b 
for  streamwise  velocity,  u,  density,  p,  temperature,  T,  and  Mach  number,  M.  Under  the 
parallel  flow  approximation,  the  transverse  velocity,  v  is  assumed  to  be  zero.  The  presence 
of  a  generalized  inflection  point,  d  ( du/dy )  /T  /dy  —  0,  at  y  —  2.775  ensures  the  existence 
of  unstable  modes  whose  properties  have  been  obtained  by  Fedorov  and  Turnin.36  The 
most  unstable  mode  is  then  introduced  at  the  upstream  end  of  the  domain  and  permitted 
to  grow.  Two  different  disturbance  amplitudes  are  considered.  The  low  value, 
em  =  \(pu)'\max/ (Pu)  =  0.565  percent,  is  chosen  to  verify  the  method:  in  Ref.20  it  is  shown 
that  the  comparison  with  linear  stability  analysis  is  excellent.  The  high  value,  em  =  8.47 
percent,  is  employed  to  investigate  control  -  only  these  results  are  summarized  here. 

To  illustrate  the  growth  of  the  disturbance,  Figure  10c  exhibits  the  mass-flux  disturbance, 
(pu)f  at  a  fixed  phase  of  the  input  signal  for  the  high-forcing  case.  In  contrast  to  the 
low-forcing  case  (not  shown),  the  growth  here  is  nonlinear,  as  the  structures  distort  and 
become  elongated  in  the  transverse  direction. 

To  examine  control,  it  is  assumed  that  a  weakly  ionized  region  is  established  in  the  flow  by 
a  suitable  nonequilibrium  technique  as  in  earlier  problems.  The  area,  shown  in  Figure  11a, 
encompasses  the  region  of  significant  mass- flux  disturbances  of  Figure  10c,  and  is  chosen  to 
focus  the  electromagnetic  interaction  to  peak  growth  regions  and  to  minimize  the  impact  of 
boundary  conditions.  The  finite  streamwise  extent  also  facilitates  analysis  of  perturbation 
growth  downstream  of  the  control  region.  The  comments  made  in  Section  3  on  physical 
values  for  conductivity  apply.  Several  different  orientations  are  considered  for  the  electric 
and  magnetic  fields.  In  each,  the  R-vector  lies  in  the  x  —  y  plane  while  the  E  field, 
measured  with  the  load  factor  K  =  — E/(UB ),  points  in  the  z  direction  as  depicted 
schematically  in  Figure  lib.  When  the  magnetic  field  is  oriented  streamwise  and  the 
electric  field  is  shorted,  denoted  case  BXK0,  the  direct  effect  of  the  ponderomotive  force  is 
on  the  transverse  component  of  the  disturbance,  and  the  mean  flow  is  unaltered.  In  the 
other  cases,  the  magnetic  field  is  transverse  and  three  subcases  are  identified  based  on  the 
value  of  the  spanwise  electric  field.  ByK0  exerts  a  retarding  force  everywhere  in  the  ionized 
region.  ByK\  is  the  open  circuit  condition:  here  the  force  respectively  opposes  or  favors  the 
higher  or  lower  speed  fluid,  effectively  transferring  momentum  as  before  (Section  4).  In  the 
last  case,  denoted  ByIi2,  the  force  accelerates  the  fluid  everywhere  in  the  ionized  region. 
The  interaction  parameter  is  fixed  at  0.0001,  except  for  the  streamwise  magnetic  field  case, 
for  which  it  is  raised  to  0.01  to  account  for  the  fact  that  the  pertinent  velocity  here  is  the 
much  smaller  vertical  disturbance  velocity. 
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(c)  Growth  of  disturbance  -  no  control 

Figure  10:  Inviscid  Entropy  Layer  Instability  Problem  and  Results  Without  Control 
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8 


(a)  Conductivity  variation 


Figure  11:  Plasma  Environment  Imposed  to  Control  Inviscid  Entropy  Layer  Instability 
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The  effectiveness  of  control  will  be  discussed  with  statistics  accumulated  when  an 
asymptotic  state  is  achieved.  To  quantify  the  assessment  of  disturbance  growth,  Mack’s 
energy  norm,39  denoted  1 1 1 1  (not  to  be  confused  with  the  electric  field)  is  employed: 

ro o 

\\Ev\\(x)  =  \\E\\(x,y)dy,  (5) 

J  0 

where  _ 

||E||  (x,a)  =  p^  +  p^+-^y?+-^— (6) 

Note  that  the  w'  component  is  suppressed  because  of  the  2-D  nature  of  the  present 
simulations.  The  energy  norm  is  obtained  from  the  computations  by  collecting  the  needed 
statistical  data  over  several  cycles  of  the  input  perturbation,  starting  after  a 
time-asymptotic  state  is  reached. 

Figure  12a  depicts  the  variation  of  the  integrated  quantity,  \\Ey\\  versus  streamwise 
distance,  x.  The  no-MGD  case  provides  a  reference  for  the  growth  of  the  disturbance 
energy,  with  which  the  MGD  cases  may  be  compared.  Under  the  ||U||  criterion,  only  the 
streamwise  5-field  yields  reduction  in  disturbance  energy  growth.  This  result  is  consistent 
with  the  stability  analysis  of  incompressible  boundary  layers  by  Rossow.40  Although  the 
magnetic  field  in  this  case  serves  to  inhibit  only  the  vertical  velocity  disturbance,  it  is 
evident  that  this  effect  is  sufficient  to  fundamentally  alter  the  dynamics  of  disturbance 
growth.  The  most  rapid  destabilization  (relative  to  the  no-MGD  case)  is  observed  with 
retarding  forces,  ByK0  while  the  accelerator  case,  ByK2  shows  similar  behavior  but  the 
growth  rates  are  somewhat  smaller  downstream  of  the  control  region  (A"  >  250).  The 
open-circuit  condition  exhibits  the  smallest  currents  and  has  no  significant  impact  relative 
to  no-control. 

Further  insight  into  the  nature  of  the  disturbance  energy  growth  may  be  obtained  by 
exploring  the  behavior  of  individual  terms  comprising  the  energy  norm.  Figures  12b  and  c 
display  the  growth  of  the  kinetic  energy  and  temperature  fluctuation  components  of  the 
energy  norm.  Upstream  of  the  control  region,  x  <  50,  all  methods  indicate  similar  values  in 
each  component.  In  the  control  region,  the  streamwise  magnetic  flux  yields  a  decrease  in 
both  components  of  the  norm,  while  the  retarding  case  shows  increases  in  both. 
Acceleration,  ByK2,  exhibits  opposite  trends  in  disturbance  kinetic  energy  and 
temperature  components,  which  decrease  and  increase  respectively.  In  fact,  in  the  control 
region,  disturbance  temperature  component  growth  rates  are  independent  of  load  factor. 

An  examination  of  the  scales  of  the  ordinates  shows  that  the  kinetic  energy  component  is 
relatively  small  compared  to  the  other  two,  thus  explaining  the  overall  effect  observed  in 
Figure  12a. 

The  enhancement  of  the  disturbance  is  also  evident  in  the  coherent  features  of  the  flow. 
Figure  13  exhibits  the  disturbance  vorticity  for  the  uncontrolled  case  (top)  and  control 
through  acceleration,  ByK2  (bottom).  Even  though  the  disturbance  kinetic  energy  for  this 
control  technique  is  observed  to  diminish  in  Figure  12b,  it  is  obvious  that  the 
corresponding  vorticity  is  significantly  enhanced.  Whereas  the  no-control  case  shows 
vorticity  increase  in  essentially  one  sequence  of  structures,  the  control  case  shows  two  such 
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(c)  Temperature  component 
Figure  12:  Effect  of  Control  on  Energy  Norms 
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6 


No-MGD 


Figure  13:  Disturbance  Vorticity  Growth  with  and  Without  Control 


trains.  The  lower  train,  whose  genesis  is  explored  further  below,  is  particularly  strong  and 
occurs  near  the  lower  control  boundary. 

Figure  14a  exhibits  mean  velocity,  Mach  number  and  temperature  profiles  at  station  B  of 
Figure  11a.  The  mean  velocity  profiles  are  relatively  straightforward  to  interpret:  an 
accelerating  force  yields  an  enhancement,  while  the  decelerating  force  generates  a  deficit, 
ffowever,  the  Mach  number  across  the  layer  is  lower  in  both  cases,  consistent  with  the 
higher  temperatures.  Clearly,  heating  of  the  layer  through  Ohmic  dissipation  is  a  common 
feature  of  both  accelerating  and  decelerating  cases. 

The  rapid  growth  of  the  disturbance  energy  when  there  is  a  significant  differential  of  force 
and  energy  across  the  shear  layer  may  be  further  related  to  other  mean  flow  quantites. 
Figure  14b  exhibits  the  evolution  of  the  generalized  inflection  point  parameter, 
d / dy(l/Tdu/ dy)  at  the  three  stations  of  Figure  11a.  The  no-MGD  case,  and  MGD 
configurations  for  which  the  impact  is  either  stabilizing  (BXK0)  or  neutral,  ByK\  do  not 
introduce  new  inflection  points,  ffowever,  BvKq  and  ByK 2  show  several  distinct  points  of 
inflection  near  the  edges  of  the  control  domain.  The  new  localized  shear  layers  develop 
precisely  near  these  locations  and  evolve  with  independent  stability  characteristics. 

To  provide  further  insight  into  the  effect  of  force  and  heating  due  to  MGD,  the 
configurations  ByK0  and  ByK2  are  subjected  to  further  numerical  experiments,  energy 
norm  results  of  which  are  shown  in  Figure  14c.  In  the  computations  marked  No  JHT,  the 
j2/cr  term  is  omitted  from  the  E  ■  j  contribution  to  the  energy  equation,  while  the  work 
term,  u.(j  x  B)  is  retained.  It  is  immediately  apparent  that  the  disturbance  amplitude,  as 
measured  by  this  norm,  is  reduced  significantly,  and  is  approximately  the  same  as  the 
no-MGD  case.  Although  ignoring  heating  affects  the  growth  of  the  norm,  it  is  observed 
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that  this  has  relatively  minor  impact  on  the  mean  velocity  profiles  i.e.,  the  velocity 
deviation  from  the  no-MGD  case  is  qualitatively  the  same  as  observed  with  Joule  heating 
on  (Figure  14a).  The  generalized  inflection  point  parameter  further  indicates  that  the  new 
inflection  points  continue  to  be  observed  even  without  Joule  heating.  These  results  show 
that  although  the  inflection  points  arise  through  the  action  of  differential  forces,  and 
appear  to  be  necessary  for  growth,  the  heating  effect  also  plays  a  key  role  in  enhancing 
Mack’s  energy  norm,  consistent  with  the  dominant  influence  of  the  thermodynamic 
component  of  the  perturbation. 
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(a)  Mean  profiles  at  Station  B 


No-MGD 
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(b)  Generalized  inflection  point  analysis 


(c)  Effect  of  Joule  heating 


Figure  14:  Effect  of  Profile  and  Heating  on  Disturbance  Growth 
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6  Energy  management 


In  recent  efforts,  Refs., 41-43  the  MGD  energy  bypass  procedure  was  explored  in  the 
nonaxisymmetric  (or  rectangular  cross-section)  device  (ND),  shown  schematically  in 
Figure  15a,  operating  at  Mach  8.  The  concept  employs  MGD  to  extract  energy  from  the 
inlet,  thus  providing  a  source  for  on-board  functions  and  a  flow  potentially  better  tailored 
for  combustion.  Optionally,  an  MGD  accelerator  is  placed  in  the  nozzle  to  augment  thrust. 
Thermodynamic  analyses  have  established  broad  parameters  under  which  benefits  may  be 
realized  (e.g.,  Ref.44).  In  Refs.,41-43  high-fidelity  simulations  on  the  configuration  of 
Figure  15a  were  employed  to  make  several  observations.  It  was  shown  that  various  3-D 
features  such  as  swept  shock-wave  boundary  layer  interactions  (SBLI)  have  a  profound 
impact  on  operation.  The  resulting  separated  flow  and  vortical  structures  couple  to  the 
current,  electric  and  ponderomotive  force  fields  giving  rise  to  complex  3-D  features  which 
demonstrate  that  2-D  and  inviscid  analyses  are  inadequate.  The  MGD  generator  shows  the 
potential  to  efficiently  slow  down  flow  in  the  inlet,  thus  decreasing  its  length,  but  separation 
limits  the  useful  length  of  the  generator  through  constriction  of  electric  current  and 
attendant  Joule  heating.  Even  with  segmented  electrodes,  Hall  currents  cause  asymmetries 
and  result  in  yawing  moments.  These  considerations  are  now  examined  in  the  context  of  an 
alternative  axisymmetric  configuration  designed  to  accomplish  similar  objectives. 

The  axisymmetric  device  (AD)  is  shown  in  Figure  15b.  Each  component  in  the  ND  has  an 
equivalent  in  the  AD,  including  a  compression  surface,  a  cowl,  a  generator,  a  constant  area 
isolator/combustor,  a  nozzle  and  an  afterbody.  The  freestream  flow  parameters  are  set  to 
be  M  =  8,  Re  =  1.6  x  106  and  T, 'x  =  250K,  yielding  Uoo  =  2535 m/s.  The  nose-cone  angle 
is  set  to  16°,  to  obtain  a  Mach  number  reduction  from  8  in  the  freestream  to  4  in  the 
isolator/combustor  region.  The  nozzle  is  formed  by  reducing  the  centerbody  at  a  4°  angle 
in  the  aft  region  that  also  serves  as  an  accelerator.  The  structure  of  the  grid  is  shown  in 
Figure  15c  where  only  some  of  the  182,  512  points  are  shown.  The  mesh  size  is  determined 
from  previous  experience42  on  the  rectangular  configuration,  where  a  mesh  resolution  study 
has  been  presented.  Eleven  points  are  employed  in  the  azimuthal  direction,  since  even 
when  plasma-flow  control  is  invoked,  the  electromagnetic  field  described  below  ensures 
axisymmetry  of  the  fluid  dynamic  response.  However,  for  clarity,  in  some  figures  the  entire 
device  is  displayed  by  rotating  the  computed  solution  around  the  axis.  In  cases  where 
MGD-bypass  is  simulated,  an  azimuthal  magnetic  field  is  imposed,  as  shown  in  Figure  16a 
for  the  generator  section  (the  accelerator  has  a  similar  field  arrangement).  The  magnitude 
is  normalized  to  unity  at  the  surface  and  diminishes  in  the  radial  direction  as  1/r, 
consistent  with  an  axial  current.  Forward  and  return  paths  for  this  postulated  current 
must  be  provided  separately  for  each  component  through  proper  channels,  which  are  not 
modeled.  The  field  is  unperturbed  by  induced  currents  since  the  low  Ra  approximation  is 
invoked.  The  interaction  of  currents  in  the  fluid  with  the  external  circuit,  which  is  not 
simulated,  occurs  through  four  pairs  of  segmented  electrodes  mounted  as  shown  in  each 
figure.  Unlike  in  the  ND  case  of  Figure  15a,  here  the  electrodes  are  circular  rings  on  the 
inner  and  outer  surfaces.  Electrode/insulator  junctures  are  assumed  to  lie  along  grid  lines 
to  facilitate  the  imposition  of  boundary  conditions.  The  electrical  conductivity,  a  is  again 
assumed  to  arise  from  a  nonequilibrium  method  for  both  generator  and  accelerator  (see 
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(a)  Nonaxisymmetric  (rectangular)  configuration 


(b)  Axisymmetric  configuration 
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Coarsened  for  clarity 


(c)  Schematic  of  grid 


Figure  15:  MGD  Energy  Bypass  Configurations  Examined 
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(a)  Magnetic  field 


1:3:3 


Accelerator 
K=1 .2 


(c)  Potentials 

Figure  16:  Plasma  and  Electromagnetic  Environment 
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Table  1:  Cases  Computed.  Q= Interact  ion  parameter,  bill  parameter 


Case 

Heat 

Addition 

Q 

Rh 

1  (Baseline) 

No 

0 

0 

2 

Yes 

0 

0 

3 

No 

1 

0 

4 

Yes 

1 

0 

5 

Yes 

2 

0 

6 

Yes 

1 

1 

Figure  16b).  Products  of  modified  Gaussians  (MG)  are  chosen  as  described  in  Ref.41  Heat 
release  due  to  combustion  is  modeled  through  a  simple  source  term  in  the  total  energy 
equation.  The  spatial  distribution  in  a  region  surrounding  the  center  of  the  combustor 
segment,  also  depicted  in  Figure  16b,  is  obtained  in  a  manner  similar  to  that  employed  for 
electrical  conductivity.  The  dimensionless  heat  release  density,  Q*c  =  QcLref  / (pre/C/^e/)  is 
fixed  at  20  percent.  This  value  is  chosen  to  be  sufficient  to  generate  discernable  trends 
without  causing  thermal  choking  in  the  channel  even  with  the  upstream  MGD  generator 
operational. 

Several  cases,  with  different  parameters  as  summarized  in  Table  1,  are  computed  to  explore 
the  observed  interaction.  The  no-MGD  simulation,  Case  1,  serves  as  a  baseline.  The  effect 
of  heat  release  (without  MGD)  is  considered  in  Case  2  while  the  impact  of  MGD  (but 
without  heat-release)  at  interaction  parameter  of  1  is  explored  in  Case  3.  Both  effects  are 
considered  simultaneously  in  Case  4.  Interaction  parameter  variation  is  explored  with  the 
case  designated  5,  for  which  Q  =  2  and  Rh  =  0.  Finally,  Hall  effects  are  considered  in  Case 
6,  at  a  Hall  parameter  of  unity. 

The  specification  of  potentials  on  the  electrode  surfaces  depends  on  the  chosen  load  factor, 
local  velocity  and  magnetic  field.  The  load  factor  was  fixed  at  K  =  0.8  for  the  generator, 
and  1.2  for  the  accelerator.  In  contrast  to  the  simulations  reported  in  Ref.,42  where  the 
potentials  were  fixed,  in  the  present  calculations  these  were  allowed  to  vary  in  time  to 
steady  state  by  adapting  to  the  evolution  of  local  velocity  with  control.  Since  the  velocity 
variation  in  the  generator  and  accelerator  region  is  relatively  small,  the  fluctuation  in 
potentials  is  also  minor.  Steady  state  potential  values  for  Case  4  at  each  electrode  are 
shown  in  Figure  16c. 

The  results  are  first  described  in  terms  of  flow  field  structure.  Pressure  contours  are  plotted 
in  Figure  17,  where  a  significant  aspect  ratio  distortion  has  been  applied  to  the  display  for 
clarity.  The  baseline  case  shows  a  relatively  straightforward  shock  structure,  with  a 
primary  forebody  shock  reflecting  off  the  cowl  and  subsequently  undergoing  successive 
reflections  in  the  interior  of  the  domain  after  being  weakened  by  the  expansion  from  the 
cone  shoulder.  Heat  release,  Figure  17b  does  not  modify  the  initial  set  of  reflections  but 
yields  new  structures  near  the  outer  boundary  in  the  isolator /combustor  region  and  an 
increase  in  pressure.  When  the  MGD  components  are  operational,  Figure  17c,  the  effect  in 
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the  generator  is  more  obvious  than  in  the  accelerator.  The  first  shock  reflection  is  similar, 
but  near  the  nose-cone  shoulder  the  upstream  influence  of  the  shock/boundary  layer 
interaction  is  considerably  larger.  Changes  in  the  boundary  layer,  discussed  later,  include 
more  smeared  features  near  the  boundaries.  When  both  the  generator  and  the  heat  release 
are  on,  the  shock  train  persists  a  longer  distance  into  the  channel,  and  merges  with  the 
features  observed  with  heat  release  alone.  These  trends  in  the  generator  are  exaggerated 
when  the  interaction  parameter  is  increased,  Figure  17e,  and  the  features  observed  in  the 
combustor  are  less  prominent.  Note  particularly  the  increase  in  upstream  influence  near 
the  cone  shoulder.  The  Hall  effect  also  displays  similar  features  as  increase  in  interaction 
parameter,  though  there  are  differences  in  performance,  which  will  be  highlighted  later. 
Overall,  the  effect  of  MGD  is  to  increase  pressure  values  in  the  device,  with  highest  values 
being  observed  when  the  Hall  current  is  considered  in  the  analysis. 

Mach  contours  for  the  different  simulations  have  also  been  examined.45  Those  results  are 
not  reproduced  here,  but  are  consistent  with  the  pressure  variations.  Heat  release  and 
MGD  reduce  Mach  number.  The  accelerator  tends  to  compensate  for  Mach  number 
reductions  in  the  generator  and  the  combustor.  However,  smallest  exit  Mach  numbers  are 
observed  when  the  Hall  effect  is  considered  -  this  observation  is  related  to  heating  as 
discussed  later. 

Figure  18  exhibits  particle  traces  for  each  case.  In  all  simulations,  separation  is  observed  at 
the  shoulder  of  the  nose-cone,  where  it  experiences  a  combination  of  the  rapid  expansion 
coupled  subsequently  with  a  shock  reflected  from  the  cowl  leading  edge.  The  flow 
downstream  is  relatively  benign  for  the  baseline,  but  additional  separated  regions  are 
evident  in  each  of  the  other  cases.  Generally,  these  features  are  significantly  more 
pronounced  on  the  inner  surface  -  this  is  understandable  since  all  things  being  equal,  the 
magnetic  held  and  hence  the  (retarding  in  generator)  body  force  is  larger.  Combustion 
triggers  two  separation  cells  on  the  inner  body,  while  the  generator  by  itself  yields  a 
separated  region  upstream,  in  its  domain  of  operation.  Thus,  the  case  with  MGD  and 
heating,  Figure  18d,  displays  four  such  structures,  though  those  associated  with 
combustion  are  smaller.  Both  the  higher  interaction  parameter,  and  the  Hall  effect, 
accentuate  the  first  separation  region,  and  result  in  significant  blockage  of  the  how.  Since 
axisymmetry  is  enforced  in  the  simulation,  these  separated  structures  do  not  exhibit  the 
unsteadiness  anticipated  in  the  full  3-D  situation. 

Mean  prohles  of  the  velocity  and  Mach  number  through  the  howpath  are  shown  in 
Figure  19.  The  baseline  case  exhibits  the  smallest  change  in  each  quantity.  The  velocity 
falls  by  only  about  15  percent  in  the  internal  how  region,  while  the  mean  Mach  number 
reaches  about  3.5  in  the  constant  area  section,  and  rises  in  the  nozzle  and  aft  regions  to 
about  6.  By  comparing  the  appropriate  cases  relative  to  no  control,  both  heat  release  and 
MGD  generator  ( Q  =  1)  action  independently  decrease  the  Mach  number  by  about  0.8  and 
0.5  respectively.  The  combined  action  yields  reductions  of  about  1.3,  which  is  consistent 
with  the  fact  that  the  two  effects  are  spatially  independent  of  each  other  and  are  effectively 
decoupled  in  the  present  phenomenological  formulation.  Doubling  the  interaction 
parameter  reduces  the  Mach  number  by  an  additional  0.6.  The  lowest  Mach  numbers  are 
encountered  when  Hall  currents  are  on,  with  minimum  mean  values  reaching  about  1.5. 

Figure  20  depicts  static  and  total  pressures  along  the  centerline.  To  reduce  clutter,  only 
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(b)  Case  2 


(c)  Case  3 


(d)  Case  4 


(e)  Case  5 


(a)  Case  1 


(f)  Case  6 

Figure  17:  Pressure  Contours  in  Axisymmetric  Device 
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(a)  Case  1 


(b)  Case  2 


(c)  Case  3 


(d)  Case  4 


(e)  Case  5 


Figure  18:  Streamlines  for  Different  Cases  Examined 
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(a)  u-velocity 


Figure  19:  Mean  Profiles  through  Axisymmetric  Propulsion  Device 
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the  four  more  interesting  simulations  have  been  plotted.  Various  features  of  the  upstream 
flow  are  clearly  visible:  note  for  example  the  rapid  variations  across  the  shocks  arising  from 
the  centerbody.  For  the  baseline  case,  the  main  total  pressure  drop  occurs  at  the  inviscid 
shocks  while  combustion  yields  another  significant  reduction  at  X  ~  11.  With  both  MGD 
and  heat-release,  the  undulations  in  the  region  downstream  of  the  nose-cone  shoulder  are 
considerably  larger  than  when  either  effect  is  considered  alone.  This  results  from  the  larger 
separated  region,  through  which  the  centerline  (which  is  not  a  streamline)  penetrates.  Both 
the  Hall  effect  and  increased  interaction  parameter  cause  significant  additional  reduction  in 
total  pressure  while  the  impact  of  the  accelerator  is  again  generally  small.  The  baseline 
calculation  indicates  a  relatively  constant  total  pressure  near  the  centerline,  since  boundary 
effects  are  small  there.  Total  temperature  values  (not  shown)  have  also  been  examined.45 
Briefly,  these  indicate  a  sharp  rise  due  to  heat  release  and  a  more  modest  increase  in  the 
accelerator  region  downstream  (through  Joule  heating).  However,  in  the  generator  region, 
extraction  of  energy  is  efficient  enough  to  yield  a  roughly  7  percent  reduction.  In  general, 
the  effects  of  increasing  the  interaction  parameter  to  2  and  incorporating  the  Hall  effect 
(maintaining  Q  =  1)  are  similar  in  the  generator  and  combustor  region,  but  are  strikingly 
different  in  the  accelerator.  The  Hall  current  effectively  eliminates  the  velocity  increase 
anticipated  in  the  accelerator,  and  the  Mach  number  profile  is  essentially  flat.  Both 
increase  heat  deposition,  though  the  Hall  effect  causes  more  parasitic  effects  than  increase 
of  interaction  parameter. 

The  flow  held  phenomena  described  above  may  be  correlated  to  the  electromagnetic 
parameters  established  by  the  MGD  interaction.  In  most  cases,  only  Cases  4  and  6  will  be 
shown,  since  these  highlight  the  main  features  observed  in  the  simulations.  Figure  21 
exhibits  the  nondimensional  current  patterns  observed  in  the  two  cases,  with  current  paths 
on  top  and  contours  of  magnitude  below.  The  normalization  factor  for  current  is 
arefUrefBref.  physical  values  may  be  extracted  in  the  manner  indicated  in  Section  3.  The 
amplitude  in  both  cases  shows  a  degree  of  granularity,  which  depends  on  the  segmentation 
ratio.  As  expected,  in  Case  4  (Figure  21a)  the  current  hows  in  the  U  x  B  direction, 
overcoming  the  electric  held  and  performing  work  against  it.  In  the  accelerator,  the  current 
direction  is  reversed  as  the  external  circuit  performs  work  on  the  fluid.  The  current 
magnitude  is  much  higher  when  the  Hall  effect  is  considered,  Figure  21b,  especially  in  the 
accelerator  and  near  the  electrodes,  where  voltage  buildup  occurs  (see  Ref.13).  In  the 
generator,  the  current  lines  assume  a  slightly  streamwise  orientation  but  the  direction  is 
still  consistent  with  generator  operation.  In  the  accelerator  however,  the  Hall  effect  causes 
an  abnormal  circumstance,  where,  in  major  regions  the  current  is  generally  reversed  from 
its  nominal  direction,  and  in  fact  resembles  that  obtained  in  the  generator.  This  reversal 
has  implication  on  the  ponderomotive  forces,  which  are  discussed  next. 

Figure  22  exhibits  pondoromotive  force  vectors  for  the  two  cases.  Only  a  few  vectors  are 
plotted,  scaled  to  clarify  the  force  held.  Since  force  scales  as  product  of  current  and 
magnetic  held,  and  the  latter  is  fixed,  the  force  scale  and  variation  is  derived  from  the 
current  held.  Under  ideal  conditions,  a  uniform  body  force  distribution  is  desirable, 
retarding  or  accelerating  the  fluid  in  the  generator  and  accelerator  respectively.  However, 
high-gradient  low  speed  regions,  such  as  boundary  layers,  and  separated  hows  where 
velocity  reversal  occurs,  give  rise  to  nonideal  effects  that  affect  component  performance. 
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(b)  Total  pressure 

Figure  20:  Centerline  Profiles  for  Various  Cases 
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1.8 


(a)  Case  4  -  paths 


(b)  Case  6  -  paths 


(c)  Case  4  -  Magnitude 


(d)  Case  6  -  Magnitude 

Figure  21:  Current  Pattern  Established  in  Generator  and  Accelerator 
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(a)  Case  4 


(b)  Case  6 

Figure  22:  Ponderomotive  Force  Vectors 


Focusing  on  the  generator,  for  Case  4,  the  impact  of  separation  on  the  force  field  is  clearly 
evident  in  the  force  vectoring  towards  walls.  This  variation  is  substantially  higher  than  in 
the  ND  case  (see  Ref.11).  The  force  held  becomes  relatively  more  uniform  near  the  trailing 
edge  of  the  generator  suggesting  that,  again  unlike  the  ND  case,  the  length  of  the  generator 
here  is  not  limited  by  separation  for  efficient  operation.  Near  the  walls,  the  force  points  in 
the  nearly  body-normal  direction.  The  force  held  becomes  more  uniform  near  the 
downstream  end  of  the  generator,  as  the  distance  from  the  shoulder  of  the  forebody  cone 
increases.  In  the  accelerator,  the  force  is  oriented  downstream,  and  increases  in  magnitude 
with  streamwise  distance.  The  Hall  effect  results  in  considerable  degradation  of  the  force 
pattern  (Figure  22b).  In  the  generator,  the  force  held  obtains  a  downward  inclination, 
generating  a  reflecting  shock  pattern  that  ultimately  causes  separation  of  the  boundary 
layer  on  the  cowl  wall  as  well  (see  Figure  17f).  This  effect  is  similar  to  that  observed  in  the 
ND,  but  because  of  the  axisymmetric  nature  of  the  configuration,  does  not  give  rise  to 
potential  yawing  forces.  The  degradation  in  the  accelerator  is  considerably  higher.  As 
anticipated  from  the  electric  current  directions  above,  a  significant  retarding  force  exists  in 
the  accelerator.  This  is  consistent  with  the  observations,  outlined  earlier,  on  the  lack  of 
velocity  or  Mach  number  increase  in  the  accelerator. 

A  discussion  of  the  energetics  is  conducted  by  considering  the  total  term  in  the  energy 
equation,  E  ■  j,  the  work  done  by  the  force,  —j  ■  (u  x  B'j  and  Joule  heating.  These  are 
related  to  each  other  through  the  identity: 

(7) 

A  key  indicator  of  relative  efficacy  of  control  is  Joule  heating,  represented  by  the  first  term 
on  the  right  side  of  the  equation.  Since  flow  heating  is  detrimental  to  the  present  control 
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strategy,  minimal  values  of  this  term  are  desirable.  Contours  of  Joule  heating  are  shown  in 
Figure  23a  and  b  respectively  for  Case  4  and  Case  6.  By  the  earlier  scaling,  these  values 
are  normalized  by  arefU^ejB^ej.  The  ratio  formed  with  the  freestream  energy  scale, 
rhorefUref/ Lref  is  the  interaction  parameter  Q.  Since  Q  —  1  for  Cases  4  and  6,  the 
quantities  in  Figure  23  are  effectively  the  ratio  with  the  freestream  energy  scale.  Results 
with  both  cases  show  the  granularity  associated  with  higher  values  in  regions  bounded  by 
electrodes  -  this  feature  is  more  pronounced  than  in  the  ND  of  Ref.41  possibly  because  of 
the  wider  separation  of  opposing  electrodes  in  that  case.  Quantitative  values  of  heating  are 
much  higher  when  the  Hall  effect  is  present.  Regions  of  major  heating  occur  near  the  first 
electrode  in  the  generator  and  near  all  electrodes  in  the  accelerator.  In  the  first  region,  this 
heating  is  associated  with  massive  separation  just  downstream  of  the  cone  shoulder 
(Figurel8f)  and  relatively  high  values  of  current  (see  Figure  21b).  High  heating  in  the 
accelerator  boundary  layers  adjacent  to  the  electrodes  has  also  been  observed  in  the  ND  of 
Ref.41 

Line  plots  along  the  centerline  of  the  three  terms  of  Eqn  7  are  plotted  in  Figure  24  for 
cases  4,  5  and  6.  Positive  and  negative  values  of  j  ■  E  correspond  to  net  energy  extraction 
from  and  addition  to  the  flow  respectively.  For  the  two  cases  without  Hall  effects,  this  term 
is  generally  negative  in  the  generator  and  positive  in  the  accelerator.  The  work  term, 

— j  ■  (U  x  B)  also  follows  a  similar  sign  variation,  though  minor  localized  deviations  from 
this  pattern  are  evident  in  certain  regions.  Joule  heating  is  positive  everywhere,  as 
anticipated,  and  close  scrutiny  reveals  that  it  increases  with  the  interaction  parameter. 
Overall,  energy  is  extracted  in  the  generator  and  inserted  into  the  accelerator,  and  the 
interaction  is  stronger  in  the  generator  than  in  the  accelerator.  When  the  Hall  effect  is 
considered  however,  a  clear  degradation  in  performance  is  evident.  The  total  energy  term 
(Figure  24a)  shows  the  same  qualitative  behavior,  but  the  work  term,  Figure  24a,  is 
reversed  from  the  expected  behavior  in  the  accelerator.  This  correlates  with  the  previously 
noted  reversal  of  component  operation,  compounded  by  the  Hall  held  and  the 
accompanying  separation.  The  effect  of  the  energy  management  technique  on  integrated 
forces  has  also  been  analyzed  (see  Ref.45)  and  separated  into  components  associated  with 
pressure,  viscosity  and  magnetic  terms.  The  former  two  are  obtained  by  integrating  over 
the  wetted  area,  while  the  last  is  derived  from  a  volume  integration.  The  results  indicate 
that  pressure  forces  are  the  dominant  component  of  drag  (in  the  inlet)  and  thrust  (in  the 
nozzle).  Viscous  forces  are  highest  in  the  constant  area  duct,  but  are  much  smaller  than 
those  due  to  pressure.  The  body  forces  are  similar  in  magnitude  to  the  pressure  forces.  In 
the  generator,  a  reaction  drag  is  obtained  on  the  magnet,  while  thrust  is  observed  in  the 
accelerator  for  all  cases  except  when  the  Hall  effect  is  included  in  the  analysis,  consistent 
with  the  analysis  presented  earlier. 
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(b)  Case  6 

Figure  23:  Contours  of  Joule  Heating  Terms  in  the  Energy  Equation 

7  Conclusion 

The  use  of  magnetic  fields  is  explored  to  accomplish  several  different  objectives  of 
high-speed  flow  control.  The  goals  include  heat  load  mitigation,  separation  suppression, 
instability  growth  rate  modification  and  energy  management  in  propulsion  devices.  The 
interaction  of  the  plasma  environment  with  the  fluid  depends  upon  the  orientations, 
magnitudes  and  gradients  of  the  electromagnetic  parameters  relative  to  the  fluid  velocity 
and  thermodynamic  quantities.  In  general,  the  balance  between  work  done  by  the 
ponderomotive  forces  and  Ohmic  heating  is  a  crucial  factor  in  determining  the  efficiency  of 
control.  Efforts  to  leverage  fluid  phenomena  through  small  perturbations  are  more 
attractive  than  brute  force  approaches.  Thus,  in  Type  IV  shock-on-cowl-lip  interactions, 
perturbations  to  the  primary  triple  point  realize  a  milder  interaction,  but  require  the 
establishment  of  a  suitable  plasma  environment  at  a  greater  distance  from  the  body. 
Electrodeless  methods  yield  smaller  Joule  heating  than  in  arrangements  where  electrodes 
are  present,  and  thus  tend  to  be  more  successful.  Eddy  currents  arising  when  the 
boundaries  are  insulated  can  be  utilized  to  transfer  momentum  from  high-speed  regions  to 
near  wall  low-speed  regimes  -  this  technique  shows  the  potential  to  reduce  or  eliminate 
separation.  In  the  effort  to  alter  the  growth  rate  of  unstable  disturbances  in  an  entropy 
layer,  field  orientations  that  affect  the  mean  flow  give  rise  to  thermodynamic  fluctuations 
that  dominate  the  energy  norm,  regardless  of  whether  the  imposed  force  field  is  favorable 
or  adverse.  However,  when  the  held  is  oriented  to  influence  the  disturbance  quantity  alone, 
the  damping  effect  interferes  with  the  growth  process  and  results  in  smaller  growth  rates. 
For  energy  management,  energy  extraction  or  deposition  processes  necessarily  require 
current  transactions  between  the  fluid  and  the  body  through  electrodes.  New  axisymmetric 
simulations  on  air-breathing  propulsion  devices  are  described  and  compared  with  earlier 


41 


(a)  J  ■  E 


(b)  -j  •  ( UxB ) 


(c)  Joule  heating 

Figure  24:  Balance  Between  Energy  Terms  Along  Centerline 
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nonaxisymmetric  studies.  Several  aspects  are  reiterated,  including  particularly  the 
profound  effect  of  separation  on  the  performance  of  MGD  devices  and  the  superior 
performance  of  the  generator  relative  to  the  accelerator.  However,  the  side-forces  observed 
in  the  nonaxisymmetric  configuration  because  of  Hall  effects  are  eliminated  in  the 
axisymmetric  case  because  of  the  azimuthally  oriented  magnetic  held. 
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List  of  Acronyms,  Abbreviations  and  Symbols 


B 

c 

Cp 

D 

E 

Eem 

ECBMT 

F 

-1  em 

F,  G 

H 

Ha 

IL,  JL,  KL 

J 

J 

k 

K 

L 

M 

n 

P 

Pr 

q 

Q 

Qf 

R 

Rb 

Rea 

t 

T 


U 

w 

x,y,z 

X 

Z 

5 

e 

7 

P 


magnetic  induction  vector 

speed  of  light;  wave  velocity 

specific  heat  at  constant  pressure 

electric  displacement 

electric  field  intensity 

electromagnetic  energy 

Eddy  Current  Based  Momentum  Transfer 

electromagnetic  force 

flux  vectors 

flux  vector;  magnetic  field  intensity;  height 

Hartmann  number 

number  of  mesh  points 

current  density 

integrated  current 

thermal  conductivity 

electric  load  factor 

reference  length 

Mach  number 

along  normal  direction 

static  pressure 

Prandtl  number 

charge  density 

heat  flux 

flow  rate 

gas  constant  of  medium;  residual 
magnetic  force  (or  pressure)  number 
Magnetic  Reynolds  number 
time 

temperature 

Cartesian  components  of  velocity 
velocity  vector 
width  of  channel 
Cartesian  coordinates 
vector  of  conserved  variables 

sum  of  specific  internal,  kinetic  and  magnetic  energies 

small  number 

dielectric  constant 

ratio  of  specific  heats 

molecular  viscosity 

magnetic  permeability 
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UJ 

frequency 

electric  potential 

p 

density 

cr 

conductivity 

f 

shear  stress  tensor 

£,  v,  C 

transformed  coordinates 

Subscripts: 

V 

w 


x,y,z 

ref 


viscous 

wall 

components  or  derivatives  with  respect  to  x,y,z 
reference  values 


Superscripts: 

*  nondimensional  quantity 

'  vectors  in  transformed  coordinates;  derivative 
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